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Abstract 

In the context of the semiclassical theory of short periodic orbits, scar 
functions play a crucial role. These wavefunctions live in the neighbour- 
hood of the trajectories, resembling the hyperbolic structure of the phase 
space in their immediate vicinity. This property makes them extremely 
suitable for investigating chaotic eigenfunctions. On the other hand, for all 
practical purposes reductions to Poincare sections become essential. Here 
we give a detailed explanation of resonances and scar functions construc- 
tion in the Bunimovich stadium billiard and the corresponding reduction 
to the boundary. Moreover, we develop a method that takes into account 
the departure of the unstable and stable manifolds from the linear regime. 
This new feature extends the validity of the expressions. 

PACS:05.45.Mt, 03.65.Sq, 45.05.+X 
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1 Introduction 



During the past few years, research carried out on chaotic eigenfunctions has 
provided very important results. Berry and Voros ^| conjectured that, in the 
semiclassical limit, these eigenfunctions would be locally like random superpo- 
sitions of plane waves; this conjecture is supported by theorems of Shnirelman 
H and Colin de Verdiere [Q. But Heller Q found that a large number of 
highly excited eigenfunctions of the Bunimovich stadium billiard [|| have den- 
sity enhancements along the shortest periodic orbits (POs). Since then several 
studies have focused on these phenomena and led to theoretical developments 
and experimental observations such as in macroscopic billiard-shaped microwave 
cavities ||, tunnel junctions |)) and Hydrogen atoms in strong magnetic fields 

Recently, a variety of new aproaches to study the structure of chaotic eigen- 
functions have been developed |l2], [H| pL5| . One of them consists of the 
semiclassical construction of resonances with hyperbolic structure associated to 
unstable periodic orbits These resonances were studied both analytically 
and numerically. They are classically motivated constructions that take into 
account complete classical information in the neighbourhood of a PO. The so- 
called scar functions, as they have been named, can be obtained by a linear 
combination of resonances |l6[ fT?) of a periodic orbit at a given energy with 
minimum energy dispersion. Resonances span a basis and are obtained by ap- 
plying a creation operator over a vacuum state. The vacuum state (i.e. a 
resonance with no transverse excitations) is constructed with a (conveniently 
selected) Gaussian wave packet that follows a modified transverse motion along 
the chosen orbit. The modified motion, which is the result of dropping the 
pure hyperbolic one, describes a bunch of POs surrounding the chosen trajec- 
tory. For instance, the eigenvectors of the monodromy matrix evolve without 
the exponential contraction-dilation that is expected in this kind of hyperbolic 
dynamics. Instead, after one period they return to themselves (up to a minus 
sign in some cases). In the same way, the wave packet returns to itself with an 
accumulated phase that is an integral multiple of 2tt, guaranteeing its continuity. 

Billiards are among the most interesting and well-studied systems in quan- 
tum chaos. We are going to work with the Bunimovich stadium billiard. For 
billiards, the boundary is the natural Poincare section. We present a formula- 
tion of resonances and scar functions over the boundary. This reduction makes 
the calculations easier and it is a clear advantage when exploring a great num- 
ber of eigenfunctions. Of course this is not the only reason to obtain these scar 
functions expressions, since they are more than tools for investigating numer- 
ically the structure of chaotic eigenfunctions. They are the cornerstone of the 
semiclassical theory of short periodic orbits of References [[l6| |l7| . In this con- 
text it is possible to obtain all quantum information of a chaotic Hamiltonian 
system just by knowing classical information of a small number of short POs (in 
fact, the shortest ones, whose number increases at most linearly with a Heisem- 
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berg time). This is accomplished by evaluating the interaction between POs. 
In previous calculations |l7j only vacuum states of a handful of the shortest 
POs were used. For higher energies it is necessary to incorporate excitations to 
these vacuum states and also longer POs. All this can be done by using scar 
functions associated to a greater number of POs than the ones needed at low 
energy values. It is evident that working on the billiard domain will not be the 
best choice in these cases, and expressions on the boundary become essential. 
Hence, we need to develop an efficient method to evaluate them, this being one 
of the main goals of this article. 

We also go further by taking into account the non linear behaviour of the 
unstable and stable manifolds. Then, we are able to extend the validity of scar 
functions beyond the linear regime. The chosen approach is general in nature, 
and despite it is applied to the special case of the Bunimovich stadium it can 
be extended to general systems in a straight forward way. These new ingredi- 
ents do not significantly alter the construction of scar functions, which preserve 
their compact character. These new expressions provide a more comprehensive 
description of the phase space in the neighbourhood of each orbit. They are the 
other main goal of the present work. 

This paper is organised as follows, Section 2 consists of a detailed explana- 
tion of the construction of resonances on the billiard domain. Here we address 
many subtleties involved in the wave function calculations and local expressions 
are provided. Moreover, a comprehensive approach is offered for resonances 
including vacuum states and excitations in a single formulation. Though the 
construction may seem complex the great advantage is that it is a general one, 
suitable for any billiard. In Section 3 the construction is extended with expres- 
sions on the boundary. In this part explicit formulas amenable to extensive and 
high energies calculations are presented. Also a brief summary of the results 
from the previous section and their application to the construction of scar func- 
tions on the billiard boundary is developed; here, we give explicit examples. 
Finally, Section 4 is devoted to conclusions. 

2 Resonances on the billiard domain 

In this section we give a thorough explanation of the construction of the res- 
onances associated to a given trajectory 7 of length L belonging to the Buni- 
movich stadium. We define a coordinate x along the trajectory and a coordinate 
y perpendicular to it (such that y = defines the orbit). The evolution of an 
initial displacement {y,p y ) at x = into (y(x),p y (x)) at x can be obtained by 
means of the symplectic stability matrix M(x) with elements m^, 

( v( x ) \ = ( fn X \{x) m 12 {x) \ / y \ 
\ p y (x) J \ m 21 (x) m 22 (x) J \ Py J ' 
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The eigenvectors £ u and £ s of M give the unstable and stable directions associ- 
ated to the orbit. These vectors evolve as = M(x) £ u ■ After one period 
they return over themselves, (= U (L) = (-l) p e AL £„ and £ a (L) = (-1) M e~ Ai £ a , 
where // is the total number of half turns made by £„ (x) and £ s (x) during its evo- 
lution along the orbit. For a billiard this number corresponds to fx = v + N re f 
where v is the maximum number of conjugated points (in particular, for the 
stadium v is exactly the number of bounces with the circle) . N re f is the total 
number of reflections with the billiard boundary and the symmetry lines. This 
point will be addressed later when we refer to quantisation conditions. 

We are going to decompose the motion given by M(x) into a purely hyper- 
bolic and a periodic one. In order to do so we need to specify the contraction- 
dilation rate along the manifolds. This can be done by first finding one of the 
Xq points on the trajectory where the projections of £, u (xq) and £ s (^o) on y and 
p y are equal in absolute value, i.e., the unstable and stable directions would be 
symmetrical with respect to the axes. There are 1v points like these on the 
orbit. They can be found by using the fact that M Xo (the return map starting 
at x = xq) has equal diagonal elements when this condition on the eigenvectors 
is met. By means of the relation M Xo = M(xo)M{L)M(xq)~ x this condition 
can be easily implemented. 

Then, we decompose M(x) into a periodic matrix F(x) describing the evo- 
lution of the manifolds and a matrix which is responsible of the contraction- 
dilation along them. This is just an application of Floquet theorem [ p8[ : 

M(x) = F(x) exp[f(x) A K], (1) 

where K = BDB~ X with D a diagonal matrix of elements d\\ = 1 and ^22 = — 1 
and B = (£ M £ s ), i.e., the symplectic matrix transforming coordinates from the 
new axes £ u and £ s to the old ones y and p y . The real function f(x) (required to 
fulfill f(0) = and f(L) = L) can be seen as the relation between the lengths 
of (,u(x) and £ u , but the plane y-p y has no defined norm so establishing this 
relation is in general impossible. Further conditions can be imposed on f(x), 
(see Ref. However, without loss of generality, we are going to consider 

the easier choice f(x) — x for the present calculations. On the other hand 
the neutral motion given by F(x) can be obtained by its action on £ u and £ s , 
allowing us to define the set y u {x), y s (x), p u (x) and p s {x) of periodic functions: 

(j^ ) = M*) = Hx) ^ = e-f^ M(x) (2) 

and 

) = = F(x) Z, = ef^ M(x) 6- (3) 

Since F(x) is area preserving, these functions satisfy the sympletic property 
Vu(x)p s (x) - y s (x)p u (x) = £ u (x) A Zs(x) = Zu A £ s = J ( J being the unit of 
action in the y-p y plane). 
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A resonance can be essentially seen as the product of a plane wave in the x 
direction, being the semiclassical approximation for the unidimensional motion 
along the orbit, and a Gaussian wave packet in the transverse coordinate, which 
follows a dynamics without dilation-contraction along the unstable and stable 
manifolds of the trajectory. The vacuum state (semiclassically normalised to 
unity) is given by 



;rai( 



(0) _ exp{z [S(x) + y 1 T(x)/2]/h - i 0(x)/2} 

7 1 ' VT± [n(h/J) IQ^)! 2 ] 1 ^ ' 1 ' 

where T is the period of the orbit and 4>(x) is the angle swept by Q{x) while 
evolving from to x. In this expression T(x) = P(x)/Q(x), where Q(x) and 
P{x) are the components of a complex vector constructed with the stable and 
unstable manifolds. These components are obtained as: 

?g:^!:iito^- w+,6(i)=M(i) Krj-)- < 5 > 

For billiards we take a slightly modified evolution matrix (M(x)) in order to 
have a continuous Q(x); since the phase 4>(x) of Q(x) must be known in detail 
for resonances construction, continuity is a very reasonable condition to ask for. 
This can be done by means of M(x) = (—l) N ^M(x), with N(x) being the 
number of reflections while evolving from to x. In turn, this matrix can be 
constructed by using two types of matrices: 

M ^=( I [), -d M 2 (^)=(_ 2/c 1 osW J 

Mi (I) describes the evolution for a path of length I without bounces with the 
circle (the transverse momentum is measured in units of the momentum along 
the trajectory). M%{6) takes into account a bounce with the circle (6 defines 
the angle between the incoming trajectory and the radial direction). In the 
following we assume this matrix instead of the original one given for the most 
general expressions in the previous theoretical introduction. The reflections 
related phase will be included in the expressions directly. 

Taking the wave functions defined in Eq (|4|) as the vacuum state for appro- 
priate creation-anihilation operators (see the following expression results 
for a resonance with m transverse excitations: 



\Q(x)\ 



^\x lV ), (6) 



where <j){x) is the phase introduced in Eq (||), and H m (z) are the Hermite 
polynomials (H = 1, Hi = 2£, H 2 = 4£ 2 — 2, . . .). It is easy to see that ip^ is 
also a product of two functions; the solution for the motion along the orbit and 
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m excitations of a transverse Gaussian wave packet which evolves following the 
same periodic motion previously mentioned. 

Inside each family, resonances are identified by the integer number n = 
0, 1, .. ., the number of excitations along the orbit and by m = 0, 1, . . ., the 
transverse excitations. The wave number k depends on 7, n and m through the 
rule: 

S(L)/h-N b iv/2-(m+l/2)fj,n = 2irn, (7) 

which guarantees the continuity of at x = L = 0. In this expression 

S(L) — J Q p x dx is the dynamical action and /j is the topological phase; that is, 
the number of half turns made by the manifolds along the orbit. Finally, N b 
is a pure quantum phase related to the boundary conditions (see fll5|]), which 
is equal to the number of reflections satisfying Dirichlet boundary conditions 
minus the number of reflections satisfying Neumann conditions. 
Rule (Q) can be divided into the to/i even and odd cases: 

S(L)/fi - N b n/2 - fin/2 = 2n (n + m f i/2) (8) 

for m/x even, and 

S{L)/Ti- N b Tt/2- 3^7r/2 = 2vr [n + (m- l)/x/2], (9) 

for the odd case. For the desymmetrized stadium billiard it results in 

S(L)/h = Lk ; N b = N s + s h N h + s v N v 

and n = N re f + v with N re f = N s + Nh + N v where N s , Nh and N v are 
the number of reflections with the stadium boundary and the horizontal and 
vertical symmetry lines respectively. Finally, v is the number of bounces with 
the quarter of the circle. The value of Sh is or 1 depending on the symmetry 
with respect to the horizontal axis being even or odd; this is equivalent for s v , 
where the vertical symmetry is considered. The set of allowed m values can 
range from or 1 (depending on mfj, being even or odd) to m < A e ff/2irh 
(where A e ff is the transverse area in which the construction is valid). For 
simplicity, we are going to consider 0(1) < m < n as a satisfactory criterion. 
The time reversal properties of the system under study can be applied to Eq 
, giving a real resonance expression. Resonances are constructed explicitely 
for the stadium billiard by assigning a semiclassical expression to each straight 
line of the orbit. The first one of these lines is the segment of 7 that begins 
in x\ = 0. Let x 2 (> x%) be the value of x where the path reaches the border 
of the stadium after the corresponding evolution. The path departing from 
X2 defines the second line and so on until x = L/2 is reached. If we define 
local coordinates (x^\y^') over each line in such a way that x^' = x inside 
the desymmetrized billiard, the expression for the jth line turns out to be [we 
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are going to use (x, y) in the following expressions, understanding that they 
represent the variables (x^\y^) of the jth line] 



2 

-(m + l/2)^(x) - F - mG ]. (10) 



JQ(*)I 

sm[ky 2 g J (x) + kx - Nb(xf)w/2 



Here Nt(xj') is defined in the same way as Nb but taking account of the bounces 
up to the point Xj (including it), the corresponding term takes into account the 
quantum phase associated to boundary conditions and this is the resonance 
counterpart of the phase considered in the quantisation conditions. Further- 
more, 

gj {x) = Re[T{ X )/2] = V l P ^^ , (11) 

where the x dependence is understood. We would like to point out that the 
factor k in the term ky 2 gj(x) in Eq ( |l0| ) is different from the factor 1/h in 
Eq (||). This is due to the fact that momenta in the expression for gj(x) are 
measured in units of hk, i.e., the momentum along the trajectory. The oscillator 
functions are defined by 

/(»>)(£) = ^- 1 /4 (2 ™ m!) -i/2 e -e 2 /2 jffm(0; 

where H m is the mth degree Hermite polynomial. The sine function in Eq ( |To|) 
was chosen in order to obtain real resonances by virtue of the time reversal 
symmetry of the system. We have taken Tx = L in Eq ([|) and considered 
J — TikR in Eq (|^), with k the wave number from the quantisation conditions 
and R the circle radius of the stadium. Finally &j(x) = N re f(x~j)ir + \(f>j + 
4>j(x) + ctj(x)], with ipj = (pj-x + <f)j-\(xj) + aj-\(xj) for j > 2 and (p\ = 0, 
follows the actual phase of Q(x). In the first term of this expression we introduce 
the jumps in the sign of Q{x) due to reflections on hard walls, that were avoided 
by using the M(x) matrix instead of the right one. In fact, this is the only place 
where this change is needed because T keeps its sign. The remaining term of 



this function, i.e., <pj + <f>j(x) + otj{x) is explained in the subsection 2.1. 

The initial phase of the resonances ipj , Eq jic|), is essential to satisfy the 
boundary conditions. We have set (f>\ — in the relations that define ipj, and 
account for the initial phase by means of two quantities, Fq and Go, which take 
values as shown in Fig. [j]. Fq accounts for the initial phase in order to obtain 
an even (odd) function when symmetry is even (odd), without considering the 
symmetry of the oscillator functions (this will be addressed by Go). In the 
case of having only one path to consider at the initial bounce, like in the first 
diagrams for both rows, the corresponding function must obey the symmetry 
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by itself. In the second case of the first row the symmetry related path is added 
but multiplied by —1 (see the final sum in this section and explanation thereby). 
Finally, in the two remaining cases (of the second row) an even function is needed 
since the function obtained by adding the symmetry related path multiplied by 
— 1 (1) is (anti)symmetric. Go takes into account the symmetry of the Hermite 
polynomials, which goes like m, the number of transverse excitations. Its value 
depend on the starting angle of the trajectory. If this angle is zero, the symmetry 
of the resonance is the one given by the F choice. In the other case, the 
symmetry of the oscillator function (i.e. Hermite polynomials) comes into play 
and in order to have the right final result wc need to consider it. This can be 
easily accomplished by the choice of the Go value according to Fig. |l|. 



Figure 1: Fq and Go values depending on the type of starting point and on 
the boundary conditions. Dashed (solid) lines correspond to the symmetrical 
(antisymmetrical) condition, i.e., s/,, s v — 0(1). Finally, arrows stand for the 
trajectories. 

The transformation of local coordinates (x^\y^ ) over the j line to (X, Y) 
coordinates that belong to the horizontal and vertical directions respectively, can 
be obtained in a simple way by means of an usual transformation. If (Xj,Yj) 
are the coordinates of the point Xj , and aj is the angle between the j line and 
the horizontal direction, (x"' — Xj,y^) = Gj(X,Y) is given by 



Finally, the family of resonances ip 7 is constructed by means of all the lines 
including symmetries (see Fig. ^ where this procedure is illustrated for one 
of the shortest periodic orbits of the stadium billiard) , this avoids unnecesary 
evaluation of reflection points over the symmetry lines, considering only those 
over the boundary: 




Gj(X,Y) = (X-X^Y-Yj) 



( 



cos(aj) 
sin(a i ) 



- sin(<x, ) 
cos(aj) 



) 
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^(X, y) = EEE /i '"i ^-[(^.0) + Gj(siX, Si Y)}. (12) 
j i=i 1=1 

where = and s ; = hi = [S iA + 5 h2 (l - 2s h )] and v t = [<5 U + 

^;,2(1 — 2s„)]. mh and m v depend on j and are specified like this: rrih = 1 or 2 if 
the line is symmetrical or not with respect to the horizontal axis respectively, and 
the same happens with the number m v for the vertical axis, although = 2 
and m v — 1 if the line goes through the origin. This is the only considered 
choice of the two possible ones (the other is = 1 and m v = 2) such that the 
right number of lines are included in the case of having symmetry with respect 
to the origin. 




Figure 2: Set of lines including symmetries employed in the construction of 
resonances associated to one of the shortest periodic orbits. Different coordinate 
sets involved in this construction are shown. 

This construction of resonances allows us to reduce calculations by exploiting 
the spatial symmetries of the system. We will explain the role of all these 
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indices and coodinate rotations. The values of the function ip 1 in the general 
coordinates (X, Y) in the domain are obtained from their corresponding ones 
in local coordinates over each line. This is the reason why in the case of lines 
having no symmetry with respect to the axis X or Y, rrih and m v adopt the value 
2 in such a way that the two symmetry related lines are included (in the case of 
showing that symmetry, and m v have value 1 and the only line considered is 
the original one). Then, the values given to s; and Si complete this process by 
reflecting the corresponding coordinate. The main point about this mechanism 
in Eq ([l^) is that when second lines are involved the boundary conditions at the 
axes have to be met also. To achieve this result they are added with positive 
or negative sign if the symmetry is even (sh or s v equal to 0) or odd (s^ or 
s v equal to 1). This can be directly verified by replacing i and I by 2 in the 
formulas for hi and vi , and also Sh and s v by their corresponding values (always 
referring to the reflected lines). To satisfy the boundary conditions when a line 
goes through the origin of (X, Y), only one reflection must be considered, this is 
the reason why this case is treated separately in our general rules of assignment 
for the rrih and m v values. 

Now that we have completed the construction of resonances on the domain 
of the billiard it is convenient to introduce the expression for the hyperbolic part 
of the Hamiltonian applied to these functions. As shown in Ref. |l^] it has a 
simple form in terms of conveniently defined creation-anihilation operators and 
also can be expressed as: 

-a{x) [V(m + l)(m + 2) ^ m+2) + 

+ ^(m-l)m^ m - 2) ], (13) 

with a(x) = s 2Mx/h f'(x) A/2 — s k f'(x) A (where f'(x) = 1, s is the sign of 
the slope of the unstable direction with respect to the y axis at the initial point 
x = xo , and M the mass of the particle) . The case m = corresponds to the 
expression given i n |lq| . This formula or its version on the boundary which is 
derived in Section |3j, is the cornerstone of the short periodic orbits theory and 
also plays a crucial role in the construction of scar functions. 

2.1 Following the phase of Q(x): a local expression for 
resonances 

We have pointed out that function &j(x) is the phase accumulated by Q(x) 
from the point x± = to x over the trajectory, followed in a continuous way 
and incorporating jumps due to rigid walls afterwards (we isolated them by 
using M(x)). In this subsection we are going to explain the calculation of the 
continuous part in detail. 



n 2 /2M 3 
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The function 4>j(x) — &rg[Qj(x)] — axg[Qj(xj)], where arg takes the argument 
of a complex number in the interval [0; 2ir), follows the phase of Q(x) after the 
bounce in Xj . We must be careful in the analysis because the phase of Qj (x) 
is divided by 2. Then, the proposed method will be to consider the phase of 
Qj(x) in an arbitrary point of the line j. The complex vectors Q(x) and P(x) 
are obtained by 

where 

\ h a 3 / V u e / \ P^yxj ) psyx-j ) j 

Here x^ means (for j > 2) that M (the stability matrix multiplied by (—l) Nb ( x 3 )) 
is evaluated after the bounce over the boundary in Xj . It can be easily verified 
that 

Qj(x) — [dj + (x — Xj)cj] exp[— X(x — Xj)] + 

+ i[bj + (x — Xj)dj] exp[A(a; — Xj)]. (14) 

The first conclusion that we can derive from this expression is that for x — > — oo, 
Q{x) is real, with sign equal to sgn(— Cj), and for x — > +oo is pure imaginary, 
with its sign given by sgn(dj). We will define the new variable x — x — xj. 
With this definition we can easily write Qj(x) — (aj + xcj) exp(— \x) + i(bj + 
xdj) exp(Ax). We are going to divide the x axis, having in mind the places 
where the real and imaginary part of Q(x) change signs. These are located at 
XRe — —o-j/cj for the real part and at xi m — —bj/dj for the imaginary one. 
It is easy to see that Q{xR e ) is pure imaginary with sign given by sgn[(bjCj — 
a,jdj)/cj] = sgn(— 1/cj) = sgn(-Cj). On the other hand, Q{xi m ) is real and 
sgn Q(xi m ) = sgn[(dja,j — bjCj)/dj] = sgn(l/d,) = sgn dj. Moreover, once the 
signs of Cj and dj are given, the order relation between XR e and xi m is specified. 
For instance, for Cj > and dj > 0, and using the fact that Oj dj — bj Cj = 1 
there results xr e = —aj/cj < —bj/dj = i/ m . 

The only idea that motivates the previous reasoning is the continuous evolu- 
tion of Q(x), without necessarily being a monotonous one. Hence, we can have 
four possible sign combinations and orderings of xr c and i/ m , as shown in Fig. 
| The same Figure also shows the evolution of Q(x) in the complex plane. 

In these figures we have labelled the three sectors of the x axis that are 
defined by xr c and xi m with —1, and 1. They are the three quadrants where 
Q(x) time evolution takes place, evolving from —1 to and then to 1. This 
labelling can be defined through the function: 

{1 if x > max(— aj/cj, —bj/dj) 
— 1 if x < min(— aj/cj, —bj/dj) 
otherwise. 
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Figure 3: Detail of the possible evolutions of the complex number Q(x). The 
curves in three out of the four quadrants depict its time evolution. xr c divides 
the x axis in domains with a positive or negative value of the real part of Q(x). 
The same applies to xi m but for the imaginary part. Numbers —1, and 1 
denote the quadrants visited (in the given order) by Q{x). 



Though rather technical, introducing rij allows us to define the resonances 
in a local form, and this is a very important feature of them. The curve rep- 
resenting Q(x) (in Fig. [3]) is only qualitative and, of course, it does not imply 
monotony in the evolution. In the case of Figures corresponding to xr £ > Xi m , 
for instance when dj < 0, the imaginary component of Q in the —1 region grows 
while x decreases in absolute value; but in some point it must begin to decrease 
in absolute value in order to pass through zero in a continuous way, following 
the signs of the quadrant. Something similar happens in the other cases for 
both asymptotic limits. The arrows over the real and imaginary axes in these 
four schemes represent Q(xr £ ) if they are over the imaginary axis and Q(xi rn ) 
if they are over the real axis. 

Lets put again the analysis in terms of the x variable over the trajectory. 
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Not all the changes of quadrant imply a jump in the phase of Q{x) given by the 
arg function. While it is true that inside of each quadrant there is no need of 
monotony, it is needed among the quadrants (this is a kind of discrete monotony 
that allows us to say that the phase of Q(x) is greater in 1 than in 0, and it 
is greater in than in -1). Lets suppose now that there has been a change of 
region (or quadrant). In the case of x > (there has been an evolution from the 
point Xj) but (j>j{x) < (the phase of Q(x) relative to the one that it had in Xj 
is lower), or x < and <j>j(x) > (there has been a backward evolution but the 
phase grows) , we are in a situation where the real positive axis has been crossed. 
Thus, we have to add or substract 2ir, respectively, in order to continuously keep 
track of the phase (because if x > there has been an evolution from Xj and the 
crossing has been in the counterclockwise or positive sense, and it is negative 
if x < 0). To summarize, if rij{x) ^ Tij(xj) and [x — Xj)4>j{x) < a phase 
otj{x) — 2ir sgn(ir — xj) must be added; otherwise atj(x) = 0. If Cj = or dj =0, 
Xj can be replaced by any other point over the j line, inside the desymmetrized 
billiard. 

Finally, we can see that <f>j(x) + a.j(x) defines the angle swept by Qj(x^') 
in a continuous way. Hence, (pj — ifj-i + (j)j—\{Xj) + a.j-\(xj) for j > 2 prevents 
the phase from jumping when changing from one line to the next one. 



3 Resonances on the billiard boundary 

In Section |^ we described the construction of resonances on the billiard domain. 
This is only suitable for calculations carried out in the low energy region, where 
these wave functions are simple. If we include longer orbits as well as higher 
energy families of the shortest ones we need to consider a reduction to a surface 
of section. Also, for calculations of matrix elements and for obtaining explicit 
expressions in terms of classical quantities when ?i — s- is essential to perform 
this reduction. As a surface of section we can take any differentiable curve £ 
with a coordinate q along it and another n orthogonal to £ at the point q (rj — 
over the curve). Consider an orbit 7, its Z crossings with this curve (taken to be 
a t Qj (j ' = I; • • • j 0j w ith angles 0j on £ and at Xj on 7) define a number / of mth 
excited wave packets as the function V'7" 1 '' is considered only over £. This is in 
fact the representation of on the section. In the case of billiards, systems 
that are bounded by rigid walls, we can take this boundary as the surface of 
section, equipped with Birkhoff cordinates, i.e., the boundary arclength q and 
the tangential momentum p at the bounce. The origin of q is at the point 
(X, Y) — (0, R) and it grows in the clockwise sense. But for Dirichlet boundary 
conditions - 7 is null to order \fn~ on this surface. Then, we can take 

^(q)^^-(x,y) (15) 
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as the representation of ip\ on the section. We are going to obtain the general 
expression for this function. In the neighbourhood of a bounce point, is 
given by the combination of 2 terms, one corresponding to the incoming path 
and the other to the outgoing one (see Fig. ||). 




Figure 4: Incoming and outgoing paths of a given trajectory at a bounce point 
on the billiard boundary. Coordinates {xj—\,yj-i) and (xj,yj) on the trajectory 
correspond to the incoming and outgoing paths; t,n are the tangential and 
normal coordinates to the boundary at the bounce point A (A is the origin of 
Xj and n). Finally q is the arclength coordinate, with value q = qj at the jth 
bounce. 

From this Figure we see that the normal n (the former general 77 coordinate 
specialized for this case) and tangential t coordinates can be related to the 
trajectory coordinates by means of simple rotations (being Xj — sin (9 j)t — 
cos(6j)n, yj = cos (6j)t + sin (6j)n and Xj-i — sin (Oj)t + cos (6j)n, Vj-i — 
- cos (6j)t + sin (6j)n). We can consider t ~ R((q - qj)/R - (q- qj) 3 /(6R 3 )) 
for the tangential coordinate, where R is the curvature at the bounce. Also, for 
the normal and tangential coordinates to the boundary at the bounce point, the 
relation n ~ —t 2 /(2R) is valid (over the boundary). For this reason we are going 
to take y ~ cos(9)t ~ cos(0)(q — qj) and x ~ sin(8)(q — qj) in expressions on the 
boundary at the lowest order in h. As stated above, we can write the general 
expression for resonances near the jth bounce as a sum of two contributions of 
the type described by Eq (||), arising from the incoming and outgoing paths, 
respectively. We can write a real expression for the normal derivative to the 
lowest order in h (we remember here Eq ( fL0|) and that the sin function vanishes 



14 



to order y/h on the boundary): 



1/4 



cos [fccos 2 9j(q — qj) 2 gj(xj) + ksm9j(q — qj) + A] (16) 



2 

where (£) is the same defined after Eq (|l^) in Section 2, but now taking 
£ = VtRiq-qj) cos 6j/\Q(xj)\ and A = kxj - N b (x+)n/2 - (m + l/2)$ J (x j ) - 
F Q - mG . 

As previously mentioned we have evaluated the expression on the boundary 
by means of the V~h order approximation in the x and y variables, this being 
a linear approximation in the boundary variable q. Note that the ordering of 
bounces and lines is as follows: q\ be the first bounce on the boundary of the 
stadium, then j = 1 is the outgoing line from it; this process extends up to 
L/2. It is important to mention that if two lines reach the bounce point, then an 
additional factor of 2 is needed in order to account for the different weight that 
these points have in "self-retracting orbits (librations)" with respect to turning 
points; in the case of rotations the overall factor can be considered part of the 
normalization condition. 

We want to obtain expressions for resonances on the boundary of the first 
quarter of the billiard. Symmetry properties must be taken into account and 
one way to do so is by considering "contributions" coming from bounces that lie 
outside the first quarter of the boundary (see Fig. The domain counterpart 
of these bounces are taken into account by Eq (|l^) of Section 2. In fact lines are 
continued outside the desymmetrised domain. But regarding boundary expres- 
sions we must explicitely introduce them. It is clear that in the semiclassical 
limit these contributions will tend to zero, but for finite Ti we must add to the 
previous expressions the same ones but evaluated at q'j = —qy 6j = —dj when 
qj ^ 0, and q'j = L/2 - qy, 6'j = -6j when q.j ^ L/4 (see Fig. |). The case "c" 
of Fig. ^| can be disregarded. 

Then, by adding up all the contributions at each bounce along the trajectory 
we arrive at the complete expression for the resonance on the boundary. This 
procedure is analogous to the one devised to construct the resonance on the 
domain by adding all contributions coming from different paths. For complete- 
ness we are going to define the hyperbolic part of the Hamiltonian acting over 
a resonance on the boundary as: 

-» k f'(i) A W(m + l)(m + 2) 5><"' +s, + (17) 

In the remaining part of this section we introduce corrections to the general 
expression for the resonance over the boundary at each bounce. They take into 
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q=0=L 




c (L/2+q.jBj) 



b (L/2-q - ep 



Figure 5: Horizontal and vertical symmetry related bounces and corresponding 
values of their Birkhoff coordinates on the stadium boundary. 

account the curvature of the manifolds that is not included in the previous linear 
approximation. Furthermore, they apply in a direct way to the already defined 
expressions for the resonances and the hyperbolic Hamiltonian. 

Eq describes the normal derivative of the resonance on a tangential 
section to the boundary at a given bounce point. But this is not enough if we 
want to make calculations with these expressions. We therefore need to include 
the way manifolds depart from the unstable and stable directions given by the 
vectors and £ s respectively. We will denote the points on the unstable (or the 
stable) manifold corresponding to the fixed point (qo,Po) by (d,p(q)) = (q — 
qo,p — po). A description of this situation can bee seen in Fig. 0. In this Figure 
p(q) is the function that defines the unstable manifold (this procedure applies 
equally well to the stable manifold). In the following we drop the subindex 
labelling the jth bounce so that making the notation easier. In order to proceed 
we need to change coordinates from the Birkhoff set to the ones defined by the 
unstable and stable directions. This task is done by the matrix B^ 1 : 




1G 



where B is the matrix with columns given by the unstable and stable vectors. 
They are evaluated at "half a bounce" and projected on the Poincare section. 
This can be obtained by means of: 

/ l/cos0 W 1 W a b \ 
\ cos(9 ) \ R/cos9 1 ) \ c d J ' 

where the first matrix projects the vectors on the section and the second matrix 
stands for subtracting half a bounce to the third one on the right. In this last 
case, columns are the unstable and stable vectors evaluated after the considered 
bounce (see previous Section for details). 

P 




(q > p ) 



q 

Figure 6: Description of point (q,p(q)) = (q — qo,p — po) on the unstable 
manifold corresponding to the fixed point at (qo,Po) on the Poincare section 

(q,p)- 

With these functions at hand it is easy to see that the direction of the 
chord that goes through the points (qo,Po) and (q,p(q)) is given by the vector 

+ fu(q)£s with f u (q) = s(q)/u(q) = a u q + f3 u q 2 + ... In order to evaluate the 
effect on the wavefunctions this description of the manifolds should be averaged, 
and we propose a heuristic method to do it. A numerically verified convenient 
choice turns out to be the following weighted average: 

/•<? n'dn' 9 1 

f u {q) ee f u (q')i^ = ^a u q + -f3 u ? + ■■■ (18) 

Of course, this can also be done for the function f s (q) which can be expanded 
in terms of q as 

f s (q) = la s q + ^f3 s q 2 + --- (19) 
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which gives a complete description of these second order effects. Then, the 
correction to the vectors t; u (q) and £s(<?) is the most direct method to include 
them in the resonances. The new vectors now depend on the boundary variable 
q and must satisfy also the normalization condition £ u (q) A £ s (q) — J- In fact, 
the normalization factor does not play any role in the wave function since, for 
instance, it cancels out in T(q). Therefore, we do not show it in the following 
formulae for the q dependent unstable and stable vectors: 

U<l) = (1 + vhWu + fu(q)U and Uq) = (i + vf.XZ. + fs(q)U (20) 

where r\ = — sgn(q u q s ) sin 2 9. The factors (1 + T]f u ) and (1 + T]f s ) are due to 
the fact that, for high values of p, projecting the unstable and stable vectors 
on the Poincare section results into a stretching of them in q direction and a 
contraction in p. Then, when modified to consider nonlinear effects by means of 
fu{q) and f s (q), a similar definition for both vectors is guaranteed by this factor. 
Finally, the way in which all this corrections enter the previous expressions is 
by new versions of Re[P(q)/Q(q)] and 4>{q) expanded in powers of q. It is 
worth to mention that these Q and P are different from those obtained in 
the previous section, the new ones are on the surface of section. Then, the 
final expressions can be obtained by taking into account that now, for instance, 
Q = {qu + fu(q) q s ) (I + 77 f u (q))+ i (q s + f s (q) (to) 0-+V fs(q)) and equivalently 
for P. Then, it is easy to see that 

~f~\ u ( p t~\\ Puqu+p s q s . g q Rq . rnl ^ / 01 -> 

9(q) - Re [^(q)) = + + 0(q ), (21) 

with g q = {a u +a s ){q^ l ~q1)-\-2{a u —OL s )\q u q s \ sin 2 9, and where the q independent 
part can be related to the expression for g(x) given previously (see Section 2). 
In that case expressions were valid in the domain but here they are given on the 
boundary (notice the different coordinate dependence). Finally, 



= 4, + \[{a s q 2 u - a u q 2 s ) + {a u - a s )\q u q s \ sin 2 9]-^-^ + 0{q 2 ). (22) 

9tt + q s 



Then, the expression for the resonance on the boundary (at each bounce) be- 
comes: 

/ld\ i/4 

<pW(q) = -k cos9 
2 

7l 



2 

— = cos{kg(q)q 2 + ksm(9)q + A] (23) 



with A the same as in Eq (16) but taking <p in place of <j> in the expression 
for $. As can be seen we have corrected the phase only since the amplitude 
corrections are not relevant. To summarise, we have modified the vectors and 
£ s in order to take into account the effect of the curvature of the manifolds in 
the resonances on the boundary. Our approach allows us to keep the expressions 
in their compact fashion. 
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3.1 Scar function: from the domain to the boundary 



Now that we have obtained a basis of resonances expressed on the billiard bound- 
ary for a given trajectory 7 and for a quantized energy E 1 (Eq (p3|)), we can 
translate this to the scar functions. They are linear combinations of resonances 
having the form |15| ]: 

^ = X> V^W£f =0 4 (24) 
3=0 

with minimum dispersion a, where 

a 2 =< (j) 7 \(H - £ 7 ) 2 |0 7 >=< ^\H 2 h \^ > . (25) 

The boundary representation only amounts to take the normal derivative in each 
term in the previous sum. Also, the second order corrections enter immediately 
through each resonance. We illustrate this by means of Figs. |?] and ||, where 
several examples of short periodic orbits and linear density plots of the Husimi 
distributions of the corresponding scar functions are displayed. In this latter 
case the solid lines that pass through each fixed point on the Poincare section 
represent the unstable and stable manifolds. Finally the relevance of the new 
formulation can be appreciated with the aid of Fig. ||, where the scar function 
in the linear approximation can be observed, in this case the one corresponding 
to orbit 2 (see Fig. fj] for reference) . 

For comparison purposes, we would like to show the scar functions on the 
domain also. The same set of orbits is used for constructing them. They can 
be seen in Figure pi 



4 Conclusions 

Resonances on periodic orbits form a convenient basis for the investigation of 
chaotic eigenfunctions. In the present paper we have extended the construction 
of resonances without excitations in the Bunimovich stadium billiard domain to 
resonances with transversal excitations and also to scar functions. A detailed 
explanation is given, providing an explicit local expression for these wave func- 
tions. We have also extended the construction of resonances and scar functions 
from the domain to the boundary of this billiard. This procedure enables dif- 
ferent sorts of applications, the most remarkable one being the possibility to 
extend calculations of eigenstates (in the context of the semiclassical theory of 
short periodic orbits ^|) well above the first low- lying ones. In principle 
this task could have been carried out on the domain also, but it would have been 
very demanding in terms of numerical effort and practically turns out to be im- 
possible. The construction of resonances and scar functions on the boundary is 
one of the main results of this paper. 
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We have also accounted for the departure of the unstable and stable man- 
ifolds from the linear regime, which is the other important main result of the 
present work. The departure is reflected in modified expressions for the res- 
onances and scar functions. The developed method is of general scope and 
not only can be applied to general systems but also to any kind of non linear 
behaviour of the unstable and stable manifolds. This is important to under- 
line, since the previous construction of resonances was related only to the linear 
motion around the trajectories, and its area of approximate validity shrinks 
as the energy grows. With the improved but still compact expressions for the 
resonances and scar functions we are able to work with a constant effective area. 

Moreover, the direct geometrical approach involved in accounting for second 
order effects due to the projection on the boundary Poincare section makes the 
expressions extremely suitable for calculations. A deep exploration of a great 
number of highly excited eigenfunctions arises as a realistic possibility due to 
the improved speed inherent to one dimensional calculations. This will allow us 
to study large sets of eigenfunctions and evaluate statistical measures. 

Finally, with these expressions several integrals made on the domain will turn 
out to be one dimensional Gaussian integrals now, simplifying the theoretical 
analyses. 
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Figure 8: Linear density and logarithmic contour plots of Husimi distributions 
for the even scar functions corresponding to the orbits displayed in Figure [7] 
(numbers below them are the same as for the orbits). A logarithmic scale of 
uniform level ratio 1/e from the maximum downwards was used in the con- 
tour plots. Different levels of gray, uniformly distributed, complete this picture. 
Solid lines passing through fixed points represent the unstable and stable man- 
ifolds. The wavenumbers are the nearest to 1000 allowed by the quantisation 
conditions. 
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Figure 9: Linear density and logarithmic contour plot of the Husimi distribu- 
tion for the even scar function corresponding to orbit 2 displayed in Figure [?]. 
The same scales and details of Fig. p] have been considered here. 
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